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We show that geometrical scaling exhibited by the px spectra measured by the CMS collaboration 
at the LHC is substantially improved if the exponent A of the saturation scale depends on pt- This 
dependence is shown to be the same as the dependence of small x exponent of F2 structure function 
in deep inelastic scattering taken at the scale pT — Q/2. 



Recently in Refs.[TJ |2] it has been shown that pt 
spectra measured by the CMS collaboration [3] at the 
LHC exhibit geometrical scaling. Geometrical scaling 
was first introduced in the context of Golec-Biernat- 
WiisthofT (GBW) model [3] of deep inelastic scatter- 
ing (DIS) in Ref.[S]. There, a reduced 7*-proton cross- 
section a~fp{Q'^,x) that in principle depends on two in- 
dependent kinematical variables: and Bjorken a;, for 
small X {i.e. x < 0.01 or so) does depend effectively only 
upon the ratio 
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Here (5sat(a^) is so called saturation momentum which is 
proportional to the transverse gluon density [1] 



xg{x, 



(2) 



with (Jo being dimensional constant that in the GBW 
model is equal to 23 mb. Since gluon density rises for 
small X like a power, saturation momentum is customar- 
ily assumed to take the following form: 



Qo {x/xq) 



(3) 



with constant A — 0.2 0.3, determined from the HERA 
the data. 

This new energy scale emerges in the models that in- 
corporate gluon saturation (SI [7], like e.g. color glass 
condensate [5J [S], although geometrical scaling itself is 
more general, and does not require saturation {i.e. gluon 
density may grow like a power for arbitrarily small x and 
the dipole-proton cross-secction needs not to go to a con- 
stant for large dipole sizes). Nevertheless, the existence 
of the saturation scale in strong interactions is by now 
well established. If so, it should also manifest itself in 
hadronic collisions, and it indeed does, as it was shown 
in Refs.[TJ |2] where a simple Ansatz, based on dimen- 
sional analysis, for the saturation momentum has been 
proposed: 
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Here W ^ y/s and Qo ^ 1 GeV sets the scale. It turns 
out that the charged particle spectra measured at the 



LHC at three incident energies: 0.9, 2.36 and 7 TeV ex- 
hibit geometrical scaling, i.e. they fall on one energy- 
independent universal curve F{t) 
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dydp\ Ql 
if plotted in terms of the scaling variable 
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with A ~ 0.27. In essence, geometrical scaling for the px 
spectra ([5| boils down to the prescription that allows to 
relate multiplicity distributions at two different energies 
W and W. If 
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then the transverse momenta at which Eq.Q holds, sat- 
isfy 



Pt =Pt 
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This formula is independent of Qo a-nd of the overall en- 
ergy scale of W or W . So the only relevant parameter 
of geometrical scaling is exponent A. 

Equations ([t]) and Q allow to rescale pt of known 
spectrum at energy W to another energy W' , and predict 
Pt spectrum at this energy, provided we know the value 
of A. In the following, transverse momentum spectra ob- 
tained that way will be referred to as reseated spectra. 

Alternatively, if we do not know A but we know spec- 
trum at W, we can find A by changing its value until 
equality Q is satisfied, i.e. until the rescaled and true 
spectra coincide within errors. 

In dipole models of DIS the quality of phenomenolog- 
ical fits is further increased provided one incorporates 
DGLAP dependence [TU] of the saturation scale ([2|. 
Furthermore, since one "measures" Qsat with a depen- 
dent probe {e.g. with virtual photon or a pt hadron), ef- 
fective saturation scale Qsat.cff acquires some dependence 
on the virtuality of the probe. These two effects may 
be conveniently accounted for by replacing A — > X{Q^). 
Indeed, for large DIS structure function F2 that is 
directly related to the saturation scale [3] behaves as: 
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FIG. 1. Dependence of Ks on from HERA (HERA data 
points [n] after Ref.[l2]). 



Power Aeff(Q) has been extracted from the HERA data 
[TT] (see e.g. recent Ref.[T2] and references therein) as 
shown in Fig. [T] In the same figure we plot an eyeball 
fit to the experimental points given by a simple function 



Aeff(Q) =0.13 + 0.1 
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An interesting question arises, whether exponent A 
that governs geometrical scaling in hadronic collisions ex- 
hibits any dependence and, if yes, whether it is similar 
to the one obtained in DIS. For px-dependent A formula 
(Is]) takes the following form: 
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A simple way to calculate approximate pt dependence 
of A is to use Eq. (|8| bin by bin in instead of an exact 
equation ( 11 ). For slowly varying \{pt) such a procedure 
should give a good first order approximation. To this end 
we choose to rescale transverse momenta of CMS multi- 
plicity spectra atW = 0.9 and 7 TeV to the reference en- 
ergy W = 2.36 TeV for some initial value of A. Next, we 
compare the rescaled spectra with the experimental data 
at W' and repeat the whole procedure until the rescaled 
and true spectra coincide ([t]). In that way we obtain 
X{pt). Since in general for W there is no data point at 
Px obtained from (|8|, we have to interpolate the refer- 
ence spectrum and its errors (in the following we neglect 
interpolation errors). The result of this interpolation is 
depicted in Fig. |3] by a grey band. In order to estimate 
the error of A(px) we add in quadrature errors of the W 
spectrum and the interpolated error of the reference W' 



FIG. 2. Dependence of A on px- Open circles correspond to 
A obtained by rescaling 0.9 TeV data to the reference energy 
of 2.36 TeV, whereas open triangles correspond to 7 TeV. Ex- 
ponent Acfi extracted from HERA depicted by full pentagons 
(HI) and stars (Zeus) is plotted in function of px = Q/2 (see 
text). Solid line corresponds to Eq.(10l taken at Q = 2px. 



spectrum at px and p'r^ respectively. Calling this effective 
error , we repeat the whole procedure solving equation 
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for A± = \ ± 5. The result is plotted in Fig[2] to- 
gether with Acff from DIS taken at the scale ~ 4p^. 
Both fit ( 10 ) and the data points are displayed. We see 



that indeed A does depend on the transverse momen- 
tum. The agreement between A extracted from the spec- 
trum rescaled from W = 0.9 TeV (blue circles) and from 
W — 7 TeV (red triangles) is a signature of geometri- 
cal scaling. In an interval from 0.5 to approximately 2.5 
GcV A rises slowly with increasing px- Interestingly, px 
dependence of A in this interval is in a surprising accor- 
dance with DIS Xes{Q) taken at pr = (9/2. For higher 
Pt data become too noisy to draw definite conclusions. 
The smooth behavior changes completely for px < 0.5 
GeV where the steep rise of A with decreasing px is seen. 
This may be a signal of an onset of a some other compo- 
nent in the production mechanism. Here, however, the 
assumption of slowly varying A breaks down and more 
numerical care is needed before quantitative conclusions 
concerning small px part can be drawn. One should also 
stress at this point that further analysis of low px ge- 
ometrical scaling requires good quality low momentum 
data. 

Final conclusion that has to be drawn from Figj2] is 
that geometrical scaling with constant A is certainly a 
good first approximation, but a mild px dependence of A 
improves substantially the quality of geometrical scaling. 
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This is depicted in FigjS] where we plot the pr spectra 
in terms of the rescaled momentum pi^ in the vicinity of 
1 GeV where the difference between constant A = 0.27 
and "running" A of Eq.(lO) is most pronounced. Black 



points and the shaded band correspond to the CMS spec- 
trum (and its interpolation) at 2.36 TeV. Blue and red 
points (connected by dashed lines) correspond to 0.9 and 
7 TeV spectra respectively, rescaled to the reference en- 
ergy of 2.36 TeV for constant A and "running" Aoff (2pt) 
of Eq.(10|. An improvement for "running" A is evident. 



For other px intervals where the difference between con- 
stant and "running" A is not large, the quality of geo- 
metrical scaling is ~ obviously - comparable. 
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FIG. 3. Multiplicity density for ^ = 2.36 TeV (black points 
and shaded band) as measured by CMS, and 0.9 TeV (blue 
points and dashed line) and 7 TeV (red points and dashed 
line) spectra rescaled to 2.36 TeV using hypothesis of geomet- 
rical scaling with constant A = 0.27 and "running" effective 
Aeff(2pT) for low Pt- Horizontal scale in GeV/c. 



Geometrical scaling in hadronic collisions is by far less 
obvious than in DIS. In DIS we have at our disposal sim- 
ple theoretical (GBW) model |4] that allows to identify 
kinematical variables relevant for geometrical scaling. In 
hadronic collisions such models exist [T5Hr7] but they 
rely on /ct factorization which has not been proven for 
soft particle production in central rapidity. Nevertheless, 
if fcx factorization is assumed, like in the recent studies of 
Refs. [ini E] , then the proportionality of multiplicity of 
produced gluons to the saturation momentum, and there- 
fore geometrical scaling - assuming local parton-hadron 
duality - can be derived in a rather straightforward way 
(see e.g [2]). Nevertheless, the exact form of the the 
scaling variable t, that in principle may depend also on 
rapidity, is to some extent a matter of educated guess. 
Luckily, for constant A some uncertainties cancel out in 
Eq.(|8]), showing that the only relevant parameter is ex- 
ponent A. 

Another notable difference between DIS and hadronic 
collisions is that in DIS we deal with totally inclusive 
cross-section, whereas in pp both hadronization and fi- 
nal state interactions play essential role. Nevertheless 
the imprint of the saturation scale Qsat is visible in the 
spectra, which means that the information on the initial 
fireball survives until final hadrons are formed. 

In this letter we have shown that the quality of ge- 
ometrical scaling improves if the exponent A becomes 
PT-dependent. We have computed this dependence by 
rescaling px spectra at 0.9 and 7 TeV to the reference en- 
ergy 2.36 TeV, however we have also checked that rescal- 
ing 0.9 and 2.36 TeV spectra to 7 TeV or 7 and 2.36 TeV 
spectra to 0.9 TeV gives qualitatively the same results. 
Not only spectra rescaled from different energies to 
the reference energy W' agree (which is the essence of 
geometrical scaling), but the px dependence of the ex- 
ponent A agrees with the dependence obtained from DIS 
Aofi(Q), at the scale Q ~ 2px- We find this last result 
remarkable, since it provides a direct link between two 
different types of reactions. 

Several points require further clarification. First of all 
new large px data of good quality will be of importance to 
test the range of applicability of geometrical scaling and 
of the discussed similarity with DIS. Also low px data, 
where hadronic A(px) deviates from the one from DIS, 
is required to see whether this deviation signals an on- 
set of a new production mechanism common for different 
energies, or whether different energies require different 
A(px) violating geometrical scaling in this region. It will 
be interesting to verify if geometrical scaling works also 
in heavy ion collisions. If so, px spectra in heavy ion 
collisions measured at different energies and at different 
centralities will allow find A dependence and impact pa- 
rameter dependence of Qsat- 
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